Free energies in the presence of electric and magnetic fields 
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We discuss different free energies for materials in static electric and magnetic fields. We explain 
what the corresponding Hamiltonians are, and describe which choice gives rise to which result for 
the free energy change, dF, in the thermodynamic identity. We also discuss which Hamiltonian is 
the most appropriate for calculations using statistical mechanics, as well as the relationship between 
the various free energies and the "Landau function", which has to be minimized to determine the 
equilibrium polarization or magnetization, and is central to Landau's theory of second order phase 
transitions. 
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I. INTRODUCTION 

One frequently needs to calculate the properties of 
materials in magnetic or electric fields using statistical 
mechanics. The procedure, of course, is to start from 
"the Hamiltonian", perform the statistical sum to get 
the partition function, the logarithm of which gives 
"the free energy" . Unfortunately, there is much confu- 
sion in the literature because different authors include 
different pieces of the energy of the electric or magnetic 
field in the Hamiltonian, generally without explaining 
precisely what they are doing. In this paper we de- 
scribe the different results that are given for the cases 
of static electric fields and static magnetic fields, show 
the relations between them, and explain which choice 
gives rise to which result for the free energy change dF 
in the thermodynamic identity. We also discuss which 
Hamiltonian, and hence free energy, is the most appro- 
priate. Finally we describe the connection between the 
various free energies and the "Landau function" which 
has to be minimized in Landau's theory of second order 
phase transitions. 

We use rationalized Gaussian units throughout, 
since this encumbers the formulae with a minimum 
number of factors for units and numerical factors. We 
also set ?i = c = 1. 



II. LAGRANGIAN AND HAMILTONIAN 

For a charged particle in an electromagnetic field the 
Lagrangian is 



1 If 
L = -mxo + gxo-A(xo)-g0(xo) + - / dx [E^ 

(1) 

where Xq and Xq are the position and velocity of the 
particle which has charge q, and the scalar potential 
(j) and vector potential A are related to the electric 
and magnetic fields by E = — V0 — dtA and B = 
V X A. This parametrization automatically ensures 
that V X E = and V • B 0. The justification for 
this Lagrangian is that the Euler-Lagrange equations of 
motion derived from it for A and xq correctly give 



the remaining two (inhomogeneous) Maxwell's equa- 
tions, with charge density and current density given 

by 



p(x) = (7(5(x - xo) 



J(x) = gxo(5(x-xo) , (2) 



as well as the equation of motion for the particle under 
the Lorentz force^. 

For statistical mechanics we need the Hamiltonian, 
and to obtain to this from the Lagrangian in Eq. 1^, we 
first have to construct the canonical momenta. There is 
no canonical momentum for (j) since L does not depend 
on dt(t>\ 4> is not a dynamical variable. For the other 
degrees of freedom, we have po = dL/dxo = rnxg + gA 
and Ha = dL/d{dtA) = V0 + dtA = -E. Hence the 
Hamiltonian is equal to 



n 



Po • xo + y" dx (IIa • dtA) - L 

1 2 
= 7^ (Po - 9A(xo)) +#(xo) 
Zm 
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(VxA)'+ni -V<P-Ua\{3) 



Applying Hamilton's equations of motion to Eq. ^ 
gives Maxwell's equations and the equation of motion 
for the particle under the Lorentz force, as required. 
Equation (^J can be simplified by noting that since 
is not a dynamical variable, we can insert the condition 



V • E = p 



(4) 



corresponding to d^L =^ —d^Ti — 0, in the Hamiltonian 
itself. Since IIa = — E, the last term in Eq. Q can be 
rearranged as 



dx V0-E 



dx0V-E 



dyicjyp = -g(/)(xo) , 



in which we integrated by parts in the first step, used 
Eq. Q in the second step, and used Eq. jSJl in the last 
step. Hence this term cancels the +q(j)(xo) term in the 
Hamiltonian, Eq. The Hamiltonian can therefore 



2 



be simplified to 

^= 2^(Po-gA(xo))'-7S-B + i 



, (5) 

in which the condition V-E = p must be imposed, and 
we have included the magnetic coupling to the parti- 
cle's spin s with 7 being the gyromagnetic ratio. 

In the next two sections we consider the problem 
of a many-particle system of charged particles in the 
presence of an external field produced by some external 
apparatus. An important task will be to separate out 
which pieces of the Hamiltonian should be associated 
with the system and which with the apparatus. 



III. ELECTROSTATICS 

In this section we neglect the effects of the magnetic 
field and set 



A = (soB = 0), VxE = 0, E: 



(6) 



where cj) is independent of time. We consider a situation 
in which there is a dielectric which has bound charge 
density pi, giving rise to an electric field Eb, as well as 
free charges Pext, entirely outside the dielectric, which 
give rise to an external field Eoxt- We have 



p = Pb + Pcxt : 



E 



E. 



(7) 



and the relation between charge density and corre- 
sponding electric field is given by Gauss' law 



V-E = p, 



V-Eft 



Pb : 



V • Eoxt = Pcxt • (8) 



To obtain the Hamiltonian of the dielectric, we re- 
move (1/2) /dxEp^j from Eq. (jsj, since this is the 
work done to set up the free charges, and should be 
considered part of the energy of the external appara- 
tus, not the dielectric. Hence the Hamiltonian of the 
dielectric is given by 



Hi + 7^2 , 



(9) 
(10) 



where Tio is the part of the Hamiltonian not explicitly 
dependent on the electric field, Hi is given by 



and 



Hi = / dx Eb • Eoxt , 



H2 = ^ I dxEl 



(11) 



(12) 



Note that there is an implicit dependence of Hq on the 
electric field because the interactions in Ho depend on 



the locations of the particles, and these affect the field 
through Gauss' law Eq. ((Sj). 

It is useful to express Eqs. Ijlll) and 1)12(1 in a different 
way. Writing Eoxt = — V0oxt and integrating by parts, 
Eq. ((TT|l can be written 



Hi = 



dx 



(13) 



showing that Hi represents the interaction between 
charges in the dielectric pb and the external electric 
potential (j)cxt- The integral in Eq. (|13|l is only over 
the volume V of the sample, and in the rest of the 
paper we will find it convenient to indicate integrals 
restricted to the volume of the sample by a subscript 
V. Integrals with no such indicator will be over all 
space. Similarly Eq. H12|) can be written 



H2 



I I dx 



Pfc(xi)pb(x2) 



V 47r|xi 



X2 



(14) 

showing that H2 represents the self-interaction of 
charges in the dielectric. In the second expression in 
Eq. ((TUl we used Coulomb's law. 

It is conventional to describe the response of a dielec- 
tric to an external field by its polarization P, satisfying 
V ■ P = ~pb. The "electric displacement field" is de- 
fined by D = E -I- P, so that V • D = poxt- Now any 
vector field (with appropriate boundary conditions at 
infinity) can be written^ as the sum of a longitudi- 
nal part (which has zero curl) and a transverse part 
(which has zero divergence) . Comparing the equations 
for V ■ P and V • D with with Eq. |(SJ), we have 



-E, 



Eo 



(15) 



where P^ and T)^ are the longitudinal parts of P and 
D respectively. However, P ^ — E;, (and D 7^ Eoxt) be- 
cause P also has a transverse component P"^, whereas 
the electric field is purely longitudinal since V x E = 
in an electrostatic situation. Note that even with the 
physical condition that P must vanish outside the ma- 
terial, P"^ is ambiguous, although it does not result 
in any (bound) charges and hence has no measurable 
consequences. If the polarization is a local effect, with 
dipolc moments being induced in response to the local 
electric field, one can impose the additional condition 
P = /(E) (in the simplest linear function) to 

determine the polarization uniquely. However, when 
polarization is a non-local process, e.g. in a conduc- 
tor, it is genuinely ambiguous. For example, for an 
infinitely long cylindrical conductor of radius a, one 
can add an azimuthal polarization depending only on 
r [i.e. Pe — const. g{r) for < r < a and P = for 
r > a], without changing the electric field. 

Now the integral over all space of the product of a 
longitudinal field and a transverse field is zero, which 
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is easily proved by writing the longitudinal field as the 
gradient of a scalar potential and integrating by parts. 
Hence, since P = outside the sample, 



dxP • Eext = / dxP • Eext = / dxP^ • Eext 



rfxEb-Eext. (16) 

Consequently Hi in Eq, (fTTTl can be reexpressed as 



dxP • Ep 



(17) 



Hence, putting all this together the Hamiltonian of 
the dielectric is given by 



id 



Ho 



dx P • Ep 



V 



1 I dx, / dx, ^(flM^ . (18) 
2 7y Jv 47r|xi-X2| 



The free energy is then given by 

Fdioi - -fcsTlnTre-''^^-' 



(19) 



where the trace is over the states of the system. 

In equilibrium, the free energy is minimized with re- 
spect to pb- Let us now change the external field by 
some small amount (5Eext- The bound charges will re- 
spond but this causes no change in Fdioi to first or- 
der because -Fdioi is at a minimum with respect to pb- 
Hence the change in -Fdici comes entirely from the de- 
pendence of Hi on Ecxt, (see Eq. (fT7|) ') which gives^ 



- / dxP-(5Eext: 

V 



(20) 



which is a thermodynamic definition of the polariza- 
tion. The integral is, of course, only over the sample, 
not all space, because P is only non-zero in the sam- 
ple. If instead one moves the dielectric with Epxt fixed, 
Eq. computed in the reference frame of the dielec- 
tric also gives the work done in moving the dielectrio^. 
In situations where P (x) can not be defined unambigu- 
ously, it may be more convenient to replace P by — Eh 
and integrate over all space. 

Frequently, e.g. Pippard^, a different result for the 
change in free energy is given, but this corresponds 
to a different free energy which omits the piece Hi 
of the Hamiltonian, the interaction between the dipole 
moment of the system and the external field. Denoting 
the modified Hamiltonian by H^ici have 



H'a 



diel 



Ho + l I dxE: 



(21) 



which consists of all terms that do not explicitly in- 



volve the external electric field. Now Hdu 



H' 



Jy dxP- Ecxt , and the relationship between fdioi and 



i^^icj is the same^, 



^dicl — ^dicl ^ / dxP ■ Ecxt , 

Jv 



(22) 



i.e. a standard Legendre transformation, and so, using 
Eq. (j231l, we haveSil 



^^^di, 



dx Enxt, • ^P . 



(23) 



Thus -Fjici (which depends upon the position of the 
charges) should be regarded as a function of the po- 
larization P, and its derivative with respect to P gives 
the external electric field. 

It is also interesting to consider the free energy cor- 
responding to including just the Hamiltonian Hq. Sim- 
ilar arguments to those given above yield 



6Fo 



dxE-JP. 



(24) 



The last free energy that we consider in this section 
is the combined free energy of the external apparatus 
and dielectric. From Eq. (|SJ| and Eq. @, this is given 

by 

i^comb - Fo + i y dxE2 = Fdici + i y dxE^^t . (25) 
Using Eq. H2U|) . a small change in the fields gives 



comb 



dx (Eext - P) • <5Eext 
dx (D - P) • ^Eext - 



dx E • (5Ep 



dxE-(5D. 



(26) 



did 



At the first step, the integral of P • SEcxt is extended 
to cover all space. At the second and last steps, we use 
Eq. (|15|1 and the fact that SEcxt and E are longitudinal. 
Equation H26|l is the total work done on the free charges 
in setting up the field in the presence of the dielectric 
and is the form usually presented in the E&M books 
such as Jackson^ or Grifhths^. 

Eqs. 1(201), (123, ^ and ^ describe different ther- 
modynamic derivatives in the presence of an electric 
field, each of which involves a different free energy de- 
pending upon which pieces of the electric field energy 
(1/2) / dx (El + 2Eb • Ecxt + E^^t) are included. These 
are summarized in Table When the definition of P, 
^P and D are ambiguous, one can replace them with 
Efc, (5Eb and Ecxt respectively in the final column of 
Tabled with the integrals over all space. This is be- 
cause of Eq. H15() and the fact that (5Ecxt, Ecxt and E 
are longitudinal. Of the free energies in Tabled i^did 
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I— 1 Q TTl 1 H" 1 Q 

JlcLlIlll LUlildll 


Free energy change 


Hdici = Ho + 5 / dx (Ef + 2Ei, • Ecxt) 
J 




- / dxP -^Eext [5J 




^J'dicl = J 


f dxEext -rJP [S*2] 


"Wcomb = 7^0 + ^ ^ dx (Eb + 2Ei, ■ Ecxt + Ecxt) 


'^-f^comb 


dxE ■ 5D [3,4] 


"Ho 




^ dxE-5P 

I' 



TABLE I: Different free energies for the electrostatic case. For each entry, the free energy is related to the given Hamiltonian 
by _F = — fcsTlnTre"'"^. We have written the electric field piece of the energy in terms of the separate contributions from 
the bound and free charges. In this way we clearly see that the different Hamiltonians and free energies just differ as to 
which pieces of the electric field energy, (1/2) / dx (E? + 2E6 • Ecxt + EL), 

are included. The most useful form is the first 
line, since it includes all pieces of the Hamiltoiuan involving the system. Furthermore, the expression for 5-Fdiei also gives 
the work done on the dielectric if it is moved in the external field. 



is the most useful since it includes everything associ- 
ated with the dielectric: (i) the internal interactions 
not involving the electric field, (ii) the coupling of the 
charges in the dielectric to an external field, and (iii) 
the Coulomb interaction between charges in the dielec- 
tric, see Eq. ((TH|) . Hence it is surprising that Eq. (pn|l 
seems to be less often quoted than Eq. H23|) . 

IV. MAGNETOSTATICS 

In this section we consider a static magnetic field 
and so set 

E = 0, B = VxA. (27) 

The "system" is a magnetic material with "bound" cur- 
rent density producing a magnetic field Bf,, and in 
addition there is a source of external magnetic field 
Bext, produced by a "free" current density Joxt which 
lies entirely outside the system. We have 

J = Jh + Jcxt, B = Bh + Boxt • (28) 

The relation between the current and corresponding 
magnetic field is given by Ampere's law, 

VxB = J, VxBb = Jfc, VxBext=Jcxt. (29) 

As in the previous section, we separate out which 
terms in the Hamiltonian are part of the material, from 
those which are part of the external apparatus. For the 
electrostatic case, the energy needed to create the ex- 
ternal field is (1/2) J dxE^^^, independent of the state 
of the system, so it was simple to subtract this off to 
get the Hamiltonian of the material. For the magne- 
tostatic case the situation is more complicated since a 
change in the magnetic field produced by the material 
generates, according to Faraday's law, a back EMF, E, 
in the coil producing the external field, and so work 



has to be done to maintain the external current and 
field. This work is done at a rate 

^ = -ydxJexfE. (30) 

Using V X Boxt = Jextj integrating by parts, and using 
V X E = — i9tB, we obtain 

SW = J dxBext -(JB. (31) 

Hence the total work done to maintain the external 
field if the total field increases from to B (due to 
changes in the system) is 

W = y^xBext-B. (32) 

This is part of the energy of the external apparatus 
(since it is that which maintains the external current) 
and so should be subtracted from the total energy when 
determining the energy of the system. In addition, 
we need to subtract the magnetic field energy due en- 
tirely to the external field, which is (1/2) /(ixB^^^ 
from Eq. ^ and -/ dxBl^^ from Eq. The net 

effect of this is to add (1/2) J dxB^^j. 

Hence, from Eq. Q, the Hamiltonian of the mag- 
netic material isSiS 

Hm.t = Ho ({p - qA}) + i y dx (B - Bext)' . (33) 

where Ho ({p — ^A}) is the Hamiltonian of the mate- 
rial in the absence of a magnetic field except that the 
momentum p of a particle is replaced by p — gA and for 
each spin s there is an additional energy —7s ■ B. (To 
avoid too cumbersome a notation we do not explicitly 
indicate the dependence of Ho on the spins.) 

Whereas the electric field in the electrostatic is not 
a dynamical variable, but is completely determined by 
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the charge density, the magnetic field is a dynamical 
variable. More precisely, to ensure V • B = 0, we write 
B — V X A and treat A as a dynamical variable in some 
gauge, e.g. V • A = which is known as the Coulomb 
gauge. The free energy of the magnetic material is 
given by 



-Fmat = -/esTlnTre 



(34) 



where the trace is over both the state of the system 
and the vector potential in the Coulomb gauge. 

It is conventional to describe the response of a mate- 
rial to a magnetic field by its magnetization M, where 
V X M = Jb and to define the auxiliary field H as 
H = B — M so that V x H = Joxt ■ Analogous to our 
discussion of polarization in electrostatics, by compar- 
ing with Eq. I|29(l we see that the transverse parts of M 
and H satisfy 



B, 



B, 



(35) 



However, the longitudinal part of M and H are not 
specified, unless M is a local function of B. When the 
magnetization is a non-local process, as in a supercon- 
ductor, the definition of the longitudinal part of M is 
genuinely ambiguousiS. 

Let the external magnetic field Boxt be changed by 
(5Boxt = V X (5Aoxt, and calculate the change in free 
energy. Although B and all parameters of the material 
will respond to (5Boxt , they cause no change in i^mat to 
first order (as in the electrostatic case). Therefore 



SF,r 



/ dx [V X (Bext - B)] • SA 

c 



dxJh ■ SAr 



dx M ■ 5Bo 



(36) 



where the second step used Ampere's law, Eq. (|29|l . and 
the third step an integration by parts with V x M = 
Jfc . Equation (|36|1 is given in de Gennesi-^ and Kittel^. 
As in the electrostatic case, this expression can also be 
used to find the work done when moving a magnetic 
material in an external magnetic field^. 

Next we discuss other choices that appear in the 
literature for the Hamiltonian and corresponding free 
energy. To a good approximation we can replace the 
trace over the vector potential by the vector potential 
field which minimizes the free energy. This amounts to 
ignoring thermal fluctuations in the vector potential, 
similar to the neglect of quantum fluctuations in the 
radiation field that is customary in elementary quan- 
tum mechanicsi^. 

Defining Fq = - k bT In Tr'e->^'^° we have 



fix (B - Bext) • (SB - 5B, 



(37) 



Since 5B — SBcxt is transverse, from Eq. (|35|l we can 
replace B — Boxt with M. Using Ea. (|36() . we obtain 



dFn = - / rfxM • SB 



(38) 



which can be used as the thermodynamic definition of 
the magnetization. 

Some authors, e.g. Callen^^ and Pippard'', take the 
Hamiltonian of the system to be 

H',nat = Ho ({p - qA}) + i y dx [B^ - Bl^,] , (39) 



rather than Eq. (|33|l . Equation (|39|l counts the energy 
to initially set up the magnetic field, (1/2) J dxB^^^^, 
to be part the external apparatus, but the work done 
by the external currents to maintain the external field 
(if some change is made to the magnetic material) is re- 
garded as energy of the material. Comparing Eq. (|39|) 
with Eq. and using Eq. (|23J), one can see that the 
free energy -F^^t corresponding to Ti'mat is given b;"^^ 



i^mat = F/nat " / M • Bext ■ (40) 



Hence, from Eq. we have^ii^ 

- / dxBext -^M. 



(41) 



The last free energy we consider in this section, is the 
combined free energy of the external apparatus and the 
magnetic material. From Eq. |(SJ) this is given by 



1 



F^comb = Fo+ 2 j dxB\ (42) 
From Eq. (|38|l a small change in the fields gives 



(Ji^comb = j + B] • (5B = J dxH SB (43) 

Equation (|43|l is the result given in the E&M books 
such as Jackson"^. 

Eqs. ijSni), (EHIl, (El) and (gSJ describe different ther- 
modynamic derivatives in the presence of an magnetic 
field, each of which involves a different free energy de- 
pending upon which pieces of the magnetic field energy 
(1/2) / dx(B2-f 2Bf,-Bext + B2^t) are included. These 
are summarized in Table ^ In situations where the 
definition of M (and hence H) is ambiguous one can 
replace M, S'M. and H with B^, SBi, and Bext respec- 
tively in the last column of TablelTll after extending the 
integrals to cover all space. This relies on Eq. H35|) and 
the fact that (5Bcxt , Bext and SB are transverse. Of the 
free energies in Tabled Fmat is the most useful since it 
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I— 1 Q "m 1 1 ^" /^Tl 1 Q Tl 

XldlllllLUlilcLll 


Free energy change 


Hmat = Ho({p - gA}) + i / dxBg 

J 




1 dxM-5Bext [5.111 

Jv 


H:,at = Ho({p - gA}) + i y dx (B^ + 2B6 ■ Bext) 




dxBext • <5M [5.7.131 


Hcomb = :^0({P - qA}) + i ^ dx (Bi + 2Bi, ■ Bext + BL) 




/ dxU SB [3] 


Ho({p-gA}) 


J\ 


dxM-5B 



TABLE IL Different free energies for tfie magnetostatic case. We have written the magnetic field piece of the energy in terms 
of the separate contributions from the bound and free currents. In this way we clearly see that the different Hamiltonians 
and free energies just differ as to which pieces of the magnetic field energy, (1/2) J dx (B^ + 2B5 -Bext +B^xt), are included. 
The first entry in the table is the most useful since it includes all the terms in the Hamiltonian associated with the system. 
Furthermore, the expression for SFmat also gives the work done on the material if it is moved in the external field. 



includes everything associated with the magnetic ma- 
terial: (i) the part of the Hamiltonian not involving 
the magnetic field (except that we replace p — > p — qA 
and add —7s • B for each spin) and (ii) the long-range 
dipole-dipole interaction between magnetic moments in 
the material. 



V. CONCLUSIONS 

The main results of the paper are summarized in Ta- 
blesHlandini For statistical mechanics applications, we 
have argued that the most appropriate Hamiltonian for 
the system (neglecting the Hamiltonian of the external 
apparatus) is the first line of these tables, i.e. 

Hdici = Wo + ^ y dx - El,,] (44) 

for the electric case and 

Hm^t = Ho ({p - qA}) + i y rfx (B - Be,t)2 . (45) 

for the magnetic case. The corresponding free energy 
changes are 

<5Fdioi = - / dxP-^Eext (46) 
Jv 

and 

SF^^t - - / dx M • ,5Be.t . (47) 
Jv 

The relation between the first line of the tables and 
the second is given, in each case, by a standard Legen- 
dre transform, see Eqs. H22|) and (|40|) . Unfortunately, 
Eqs. (|46|l and H47|l do not figure prominently in sta- 
tistical mechanics textbooks. Reference |Q is a no- 
table exception, but unfortunately is not in print. By 
contrast, for problems in electricity and magnetism we 



are generally interested in the total work done (on the 
free charges and currents), which is given by the third 
line of the tables, and fortunately, these are the results 
quoted in E&M books. 

It is interesting to note that although the final re- 
sults for the free energy changes in the electric and 
magnetic cases are exactly analogous to each other, 
there are differences in the Hamiltonians from which 
these results are derived. The first difference is that 
the field energy involves (1/2) / dx (E^ - E^^J = 
/ dx (E^/2 — Ef,-Eext) in the electric case as opposed to 
(1/2) / dx(B-Boxt)^ = (1/2) / dxBg. The extra piece 
in the magnetic case, / dx (B^^^ — Bcxt 'B), comes from 
the work done by the external sources in the magnetic 
case to maintain the external currents against the back 
EMF produced by a change in the field of the system. 
The second difference is that whereas Ho does not have 
an explicit dependence on the electric field (though it 
depends on it implicitly, since E is tied to the charge 
density through Gauss' law, Eq. (^), it does depend 
explicitly on the magnetic field, since one has to re- 
place the velocity x by p — qA where p is the canonical 
momentum. For the case of magnetism due to spins. 
Si, there is also an explicit interaction — 7X^4^* ' B in 
Ho- Classically, one can change variables from p to 
p' = p — qA, there is no spin magnetism, and the field 
energy only involves B^, so there is no magnetism^^. 
However, since in quantum mechanics there are con- 
straints on p and not p', for example that the angular 
momentum r x p must be a multiple of h, there is a 
coupling between the dipole moment of the material 
and the external field from the Ho{{p — qA}) part — 
or the spin term — of the Hamiltonian. 

It is useful to compare the equilibrium free energies 
we have discussed so far with the Landau free energy 
function which forms the basis of Landau's theory of 
second order phase transitionsi^, and which is appli- 
cable away from equilibrium. The partition function, 
required in evaluating the free energy through statisti- 
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cal mechanics, is written (for the magnetic case) as 
Z(Bext) = 5^Tr'exp[-/5H„,at] 

m 

= ^exp[-/3FL(m,Be.t)], (48) 

m 

where Tr ' is a restricted trace, carried out over states 
with a specified value of total magnetic moment m, 
and the Landau free energy Fi^ is defined through this 
equation. In Eq. H48|l. the sum is dominated by the 
value of m which minimizes i^L at fixed Bext (assumed 
independent of x here); this is the equilibrium magne- 
tization at Bcxt- Denoting this equilibrium value of m 
by m*, we haveii 

FL(m*,Bext) =F,nat(Bext). (49) 

For the case of spin magnetism, Eq. (|48|l simplifies 
further. As seen in Eq. (O, the coupling between the 
external field and the system occurs only through the 
term — m • Boxt, (m = 1^2% ^i)' the part of —m • B 
involving the external field. Therefore 

FL(m,Bext) = F(m)-m-Bext (50) 

where F(xn) does not depend on the externally applied 
magnetic field. In this case, comparing with Eq. (|40|l 
we see thaliii 

F(m) = F;,, . (51) 



A similar process can also be carried out for the 
electrostatic case, since the coupling term in Tidici is 
— J dxEb ■ Eoxt = — JdxP • Ecxt- However, for orbital 
magnetism, Boxt cannot be separated out cleanly like 
this due to the diamagnetic term, and one is limited to 
Eq. (|48|l if the diamagnetic contribution is significant. 



Finally, we observe that the 'external' parts of 
the Hamiltonian, that were discarded in going from 
Eq. lO to Eqs. ® and were (1/2) / dxE^^^ and 
— (f/2) /dxBoxt- The sign reversal between the two 
leads to the well known resulti^ that free charges move 
to minimize the energy of the electric field (the position 
of the charges in a conductor can be viewed as such a 
minimization problem, with constraints) whereas loops 
carrying free current (held constant by batteries) move 
to maximize the energy stored in the magnetic field. 



We hope that a brief discussion of the points dis- 
cussed here will eventually find their way into text- 
books on statistical mechanics so that future students 
will be less confused about the free energy of a system 
in an electric or magnetic field than were the authors 
of this paper. 
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